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I.  INTRODUCTION 


One  of  several  sources  of  ballistic  reentry- vehicle  dispersion  occurs 
during  separation  of  the  reentry  vehicle  from  the  boost  vehicle,  which  results 
in  an  angular  misalignment  between  the  vehicle  axis  of  symmetry  and  the 
velocity  vector  at  entry  (Ref.  1).  This  initial  angle  of  attack  causes  lift  as 
the  vehicle  traverses  the  atmosphere  and  can  result  in  dispersion  from  lift 
nonaveraging.  The  dispersion  depends  on  the  manner  in  which  the  lift  vector 
and  its  rate  of  rotation  in  space  (precession  rate)  are  coupled  during  the  initial 
reentry  motion  (Ref.  2).  It  has  been  shown  that  there  are  two  precession 
modes  that  comprise  the  familiar  epicyclic  motion  of  an  axially  symmetric 
missile  in  untrimmed  flight  (Ref.  3).  Limiting  cases  of  the  two  precession 
modes  correspond  with  specific  deployment  conditions  during  boost- vehicle 
separation.  The  negative  precession  mode,  which  results  from  an  initial 
angle  of  attack  in  the  absence  of  exoatmospheric  coning  motion,  results  in 
significantly  greater  dispersion  than  the  other  limiting  case  of  positive 
precession,  in  which  the  exoatmospheric  precession  cone  is  symmetric  about 
the  velocity  vector.  A simple  closed-form  solution  obtained  for  the  dispersion 
velocity  for  the  negative  precession  case  reveais  the  first-order  influence  of 
various  parameters  on  the  dispersion.  The  dispersion  is  relatively  insensitive 
to  most  vehicle  aerodynamic  properties,  with  the  exception  of  static  margin, 
and  is  independent  of  vehicle  mass  properties  except  for  the  radius  of  gyration 
in  roll.  The  dispersion  also  varies  directly  with  the  roll  rate,  in  contrast  to 
lift  nonaveraging  dispersion  from  body-fixed  asymmetries,  which  varies 
inversely  with  roll  rate.  Also  significant  is  the  observation  that  the  trajectory 
deflection  occurs  in  the  first  cycle  of  lift  vector  precession  about  the  velocity 
vector  and  lies  in  a plane  that  leads  the  plane  of  the  initial  angular  misalignment 
by  approximately  90  deg.  If  the  magnitude  and  direction  of  the  misalignment 
angle  is  known,  this  dispersion  source  can  be  partially  compensated  for  in  the 
targeting  model. 


II.  ANALYSIS 


Cross- range  dispersion  of  a spinning  missile  can  be  described  by  the 
relation  (Ref.  2). 


V(t)  = V(0)  - ± J t 

o 


L(9)  exp(i4* ) dt 


(1) 


where  V(t)  is  the  missile  transverse  velocity  in  the  cross  plane,  L(t)  is  the 
lift  force,  which  is  dependent  on  the  angle  of  attack  9,  and  4*  is  the  preces- 
sion angle  of  the  lift  vector  in  the  cross  plane.  For  a symmetric  missile 
with  constant  roll  rate,  the  angles  9 and  'l'  are  described  approximately  by 
the  moment  equations  of  motion  (Ref.  3) 


..  . 7 Mffll  • 

9 + pp4>  sin  9 - 4*  sin  9 cos  9 = — ^ — + -j — 9 


(2) 


d 

It 


M: 
_ + 


(4<  sin  9 ) + 9 4^  cos  9 - pp9  = — j-L  ^ sin  9 


p = 0 + 4*  cos  9 


(3) 

(4) 


where  M(9)  is  the  static  pitch  moment,  and  and  are  rate  damping 
moments.  We  consider  first  the  case  of  a static  moment  only,  and  assume 
for  M(9)  the  form 


M(9) 

I 


2 

-u> 


sin  9 


(5) 


jjt 

Lateral  velocities  are  small  enough  that  the  pitch  angle  9 in  the  classical 
Euler  angle  system  is,  for  practical  purposes,  the  total  angle  of  attack. 


T 
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where  uu  is  the  vehicle  undamped  natural  pitch  frequency.  In  absence  of  a 
yaw  damping  moment  , Eq.  (?)  can  be  integrated  to  yield 

• 2 

4*  sin  8 + pp  cos  9 = const.  (6) 

Equation  (6)  is  an  expression  of  constancy  of  angular  momentum  about 
the  velocity  vector  because  there  is  no  moment  component  in  this  direction  in 
the  classical  Euler  angle  coordinates.  For  a rolling  reentry  vehicle  with 
initial  angular  misalignment  8q,  the  constant  in  Eq.  (6)  is  either  pp/cos  8q 
or  pp  cos  8q  for  the  limiting  cases  of  postivie  (4*+)  or  negative  ( tp  ) preces- 
sion, respectively,  and  Eq.  (6)  can  be  written 

, O 4 

i sin  0 = pp[(cos  8 )*  - cos  8]  (7) 

T , - O 

where  the  upper  sign  in  the  exponent  corresponds  to  positive  precession. 

The  two  limiting  cases  of  entry  angular  motion  as  8 -»9o  in  Eq.  (7)  are  shown 
in  Figs.  1(c)  and  1(d)  for  the  positive  and  negative  exponent,  respectively. 
Positive  precession  about  the  velocity  vector  (clockwise  looking  forward,  in 
the  same  direction  as  the  roll  rate)  is  initially  at  the  exoatmospheric  coning 
rate  ^ + = pp/cos  8q.  Negative  precession  (counter  clockwise)  is  initially 

zero  and  results  from  the  gyroscopic  effect  of  the  aerodynamic  static  moment 
acting  on  the  spinning  vehicle.  The  quasi- steady  approximation  0 = 8 = 0 in 
Eq.  (2)  is  a good  approximation  to  the  initial  angle-of-attack  convergence  of 
the  circular  coning  motions,  which,  with  Eq.  (5),  yields 

*+,-  =l£h  ( 1 - '/l  + 5 cos  8 ) (8) 

2 2 . 

where  p = pp/2  and  5 = u?  /p  • If  we  substitute  ^ from  Eq.  (8)  in  Eq.  (7), 
r r • , - 

we  obtain  the  relation 
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(9) 


1 2 t — 2 

2 cos  9 (cos  9 ) T = i + cos  9 ± >1+5  cos  9 sin  9 

o 

which  describes  the  quasi-steady  angle-of-attack  convergence  as  a function 
of  the  static  moment. 

2 

With  the  small  angle  approximations  sin  9 « 9 , cos  9 *=>  1 - 9 /2, 

Eqs.  (8)  and  (9)  reduce  to  the  simple  expressions  for  angle-of-attack  conver- 
gence and  precession  rate,  respectively,  of  circular  coning  motion  in  absence 
of  damping  (Ref.  3) 


J-  = (1  + 5)' 1/4  (10) 

' o 

4/  =P,(i  ±y  i + 5)  (ii) 

Equations  (10)  and  (11)  are  compared  with  (8)  and  (9)  in  Fig.  2 and  are 
found  to  be  good  approximations,  even  for  very  large  angles  of  attack.  With 
these  approximations,  we  can  express  the  net  transverse  velocity  from  Eq.  (1) 
in  terms  of  5 as  the  independent  variable  with  the  use  of  the  exponential 
atmosphere  approximation 


p(z)  = pQexp(z/H) 


(12) 


where  z is  measured  downward  from  some  initial  reference  altitude,  where 

the  density  is  p , and  H is  a constant  scale  height.  Assuming  that  the  pitch 

o 

frequency  varies  only  with  density  over  the  altitude  range  of  interest,  and 
assuming  a constant  flight  path  angle,  we  obtain  from  the  definition  of  ? 

d£  = u_sin^  dt  (13) 
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Fig.  2.  Comparison  of  Small-Angle  Approximations  for  6 and  0 
with  Exact,  Large -Angle  Values  for  = 60  Deg 
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The  lift  force  per  unit  mass  L(0)/m  in  Eq.  (1)  can  also  be  written  in 
terms  of  % according  to 


C qS9 

U0)  _ -4 


m m 


(14) 


where  f = I/mx^  is  a characteristic  length.  Substitution  of  Eqs.  (10),  (13), 
and  (14)  in  Eq.  (1)  yields,  for  the  transverse  velocity  caused  by  initial  angular 
mi  salignment 


v(5) 


if  Hp^e 


u sin  y 


2 f 

v •/ 


1 + 5)'1/4expfiA4/(C)]d5 


(15) 


where  V(0)  = 0,  and  A4M§)  is  obtained  on  integrating  from  Eq.  (11)  with 

' t “ 

the  relation  Eq.  (13).  It  is  expedient  to  make  a further  change  of  variable 
defined  by 


x = s/l  + 5 


(16) 


which  yields,  for  V(x) 


V(x) 


2if  Hp^Q 
u sin  y 


r x 

nTx 


expfi  A 4j(x)]  dx 


(17) 


The  precession  rate,  Eq.  (11),  expressed  as  a derivative  with  respect 
to  the  independent  variable  x with  the  use  of  Eqs.  (13)  and  (16)  can  be  written 


d4<  , 

' > 
dx 


2Hpr 
u sin  y 


(18) 
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The  negative  precession  rate  can  be  integrated  between  the  limits 
1 to  x to  yield 


A 4*  (x)  = -K[ln2  - 1 + x - *n(l  - x)]  (19) 

and  the  positive  precession  rate  (because  of  the  singularity  at  x = 1)  can  be 
integrated  between  the  limits  1 + e to  x to  yield 

A4*+(x)  = >^+(x)  - i|j+(1  + e)  = K[x  - 1 - fn(x  - 1)  - e - (ne]  (20) 


whe  re 


2Hpr 

u sin  y 


(21) 


The  complex  transverse  velocity,  which  determines  the  magnitude  and 
direction  of  the  trajectory  deflection,  is  determined  from  the  average  value 
of  the  integral  in  Eq . (17)  with  appropriate  integration  limits  depending  on  the 
precession  mode.  The  integral  is  a function  only  of  the  constant  K in  the 
precession  angle,  Eq.  (19)  or  (20),  in  addition  to  the  independent  variable  x. 
As  shown  in  Appendix  A,  the  nondimensional  integral 


,jf  = - iK \/~x  expfiAij' (x)]dx 


(22) 


is  only  weakly  dependent  on  K over  a wide  ange  of  K of  practical  interest. 
The  real  and  imaginary  components  of  ,9,  obtained  numerically,  for  various 
values  of  K with  both  positive  and  negative  precession  modes  are  shown  in 


Equation  (17)  is  integrated  between  the  same  limits  for  the  positive  preces- 
sion angle. 
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Figs.  3 and  4.  The  negative  precession  case  results  in  appreciably  greater 
dispersion,  and  the  integral  defined  by  Eq.  (22)  with  the  negative  precession 
angle,  Eq.  (19),  has  an  essentially  constant  average  value  .?=  - 2.0  over  a 
wide  range  of  K.  If  we  substitute  this  value  for  Eq.  (22)  in  Eq.  (17),  with  the 
definition  of  K,  Eq.  (21),  the  net  transverse  velocity  for  negative  precession 
is  described  by  the  simple  expression 


V 


= - 2fp  0 
r o 


^o 

"“st 


(23) 


The  latter  form  of  Eq.  (23)  indicates  that,  to  the  first  order,  the 
trajectory  deflection  is  a function  only  of  the  initial  angle  of  attack,  roll  rate, 
radius  of  gyration  in  roll,  and  vehicle  static  margin.  The  net  transverse 
velocity  increment  is  in  the  negative  real  direction,  which  leads  by  90  deg  the 
plane  of  the  initial  angle  of  attack  for  the  zero  coning  motion  entry  condition 
that  results  in  negative  precession.  The  trajectory  deflection  is  directly 
proportional  to  roll  rate  in  Eq.  (23),  in  contrast  to  roll- trim  dispersion  of  a 
spinning  missile,  in  which  the  dispersion  is  inversely  proportional  to  roll 
rate. 

The  simple  result,  Eq.  (23),  is  compared  in  Fig.  5 with  an  exact 
solution  for  the  real  component  v of  dispersion  velocity  during  negative  pre- 
cession, obtained  from  a numerical  integration  of  the  equations  of  motion, 
Eqs.  (2)  through  (3),  for  the  vehicle  and  trajectory  parameters  listed  below. 


h = 300  kft 
o 

I 

X 

= 1.  15  slug-ft 

u =23  kft/sec 
o 

I 

= 15  slug-ft^ 

Y = 30  deg 

x . 

St 

= 0.  2 ft 

m = 8.  703  slugs 

CA 

= 0.  1 

d = 1.  5 ft 

CXT 

= 1.9 

ty 

N0 

S = 1.767  f r 

cN(e) 

= C sin  e 

N0 

p = 2 rps 

0o 

= 20  deg 
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Fig.  3.  Real  and  Imaginary  Components  of  the  Integral  .9  with 
Positive  Precession 


An  exact  solution  for  the  imaginary  velocity  component  w is  shown  in 
Fig.  6.  A standard  atmosphere  density  distribution  was  utilized  in  the 
numerical  evaluation,  which,  in  the  altitude  range  of  220  to  300  kft,  where  the 
trajectory  deflection  occurs,  is  represented  by  an  exponental  scale  height  of 
approximately  H *»  20  kft.  This  gives  a value  of  1.98  for  K.  With  the  pre- 
ceding vehicle  and  trajectory  parameters , the  transverse  velocity  from 
Eq.  (23)  is  -2.90  ft/sec,  which  compares  with  an  average  value  of  -2.60  ft/ 
sec  for  the  velocity  component  v in  Fig.  5.  The  small  discrepancy  is 
attributed  to  the  small  angle  approximations  used  in  the  derivation  of  Eq.  (23). 
We  can  account  for  larger  angles  of  attack  by  including  in  the  integrand  of 
Eq.  (22)  the  factor  C^g/C^g,  which  is  a function  of  the  angle  of  attack  de- 
fined approximately  by 


(24) 


If  we  express  the  angle  of  attack  in  terms  of  the  independent  variable 
x using  Eqs.  (10)  and  (16),  then  the  correction  to  the  integrand  of  Eq.  (22) 

£or  CL9/CN9  i8 


(25) 


Equation  (22)  was  integrated  with  this  correction  for  K = 1.68  cor- 
responding to  the  aforementioned  vehicle  and  trajectory  parameters,  and  the 
average  value  was  found  to  be  -1.78  compared  with  -2.0  without  the  correc- 
tion. This  gives  a net  transverse  velocity  of  -2.  59  ft/sec,  which  compares 
more  favorably  with  the  exact  result  of  Fig.  5 . With  the  small  angle  approxi- 
mation of  unity  for  the  cosine  term  in  Eq.  (23),  a slightly  better  approximation 
to  the  dispersion  velocity,  Eq.  (23),  is  given  by 
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Fig.  5.  Comparison  of  Approximation  for  Dispersion  Velocity  with 

Exact  Value  Obtained  from  Numerical  Integration  of  Equations 
of  Motion 
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which  yields  -2.75  ft/sec  for  the  dispersion  velocity  of  Fig.  5. 

The  cross- range  impact  dispersion  CR  is  approximately  equal  to  the 
product  of  the  trajectory  deflection  angle  with  the  path  length  from  the  point 
of  trajectory  deflection  to  impact.  For  a straight-line  trajectory  with  con- 
stant-path angle  y,  and  with  the  net  transverse  velocity  given  by  Eq.  (23), 
the  cross-range  dispersion  is 


CR  « 


V-hl 
u sin  y 


1 P0  h4 

x o 1 

mx  u sin  y 
st  1 


(27) 


where  hj  is  the  altitude  at  which  the  net  trajectory  deflection  occurs.  We 
can  obtain  an  approximate  expression  for  h^  from  the  observation  that  the 
net  transverse  velocity  increment  is  reached  when  the  lift  vector  rotates 
through  the  first  quarter  cycle  for  the  negative  precession  case.  From  the 
definition  of  the  independent  variable  x,  Eq.  (16),  and  the  negative  preces- 
sion angle,  A4>_(x),  Eq.  (19),  the  density  at  this  altitude  is 


P1 


u2Sx 


st 


(28) 


whe  re 


t n 


x i <**  • 


2K 


- 1 


(29) 
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Again,  using  the  exponential  atmosphere  approximation,  we  obtain,  for 
hj,  in  terms  of  an  arbitrary  entry  altitude  h^,  where  the  density  is 


h . = h - Hf  n 
1 o 


p C u Sx 

' O Ng  st. 


(30) 


which  gives,  for  CR,  the  parametric  relation 


CR  « 


e ik 

o 


mllx 


st 


h - Hfn 
o 


„2/  2 \ . 2 

K ^Xj  - III  sin  y 


4p  H2 C..  Sx 

'O  Ng  St 


(31) 


with  Xj  defined  by  Eq.  (29).  The  altitude  hj  is  shown  in  Table  1 for  a range 
of  K values  and  is  relatively  insensitive  to  K.  Also  shown  is  the  cross-range 
dispersion  per  unit  angle  of  attack  0 . Hence,  the  parametric  dependence  of 
the  cross- range  dispersion  is  described  to  the  first  order  by  t’  e multiplier 
of  hj  in  Eq.  (27). 


Table  1. 

. Altitude  of 

Trajectory  Deflection 

K 

hj,  kft 

CR/0  , ft/rad 
o 

0.5 

264 

56.9 

1.0 

255 

1 10 

2.0 

245 

212 

4.0 

2 34 

403 

8.0 

221 

763 
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Aerodynamic  damping  was  excluded  from  the  equations  of  motion, 

Eqs.  (2)  and  (3),  in  the  foregoing  discussion.  In  Appendix  B,  it  is  shown 
that  damping  is  so  small  at  the  high  altitudes  where  the  trajectory  deflection 
occurs  that  it  has  a negligible  effect  on  the  dispersion. 
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III.  SUMMARY  AND  CONCLUSIONS 


A simple  relation  has  been  derived  for  the  prediction  of  cross-range 
dispersion  of  a ballistic  reentry  vehicle  that  is  caused  by  an  initial  angular 
misalignment  between  the  vehicle  axis  of  symmetry  and  the  velocity  vector  at 
entry.  The  dispersion  is  strongly  dependent  on  the  exoatmospheric  attitude 
and  motion,  which  determines  the  endoatmospheric  lift-vector  precession 
behavior.  Symmetric  coning  about  the  velocity  vector  (positive  precession) 
results  in  small  dispersion  relative  to  the  initially  nonpreces sing  vehicle  at 
angle  of  attack  (which  results  in  negative  precession).  Based  on  idealized 
first-order  aerodynamics,  dispersion  is  found  to  be  relatively  insensitive  to 
all  aerodynamic  coefficients  except  static  stability.  This  results  because 
the  initial  lift  nonaveraging  occurs  at  high  altitude  before  the  vehicle  has  ex- 
perienced appreciable  drag  deceleration  and  before  aerodynamic  damping  has 
an  appreciable  influence  on  the  angle-of-attack  convergence.  Also,  the  lift 
force  coefficient,  which  determines  the  lateral  acceleration  for  a given  angle 
of  attack,  influences  the  rate  of  angle-of-attack  convergence  for  a given 
static  margin,  and  there  is  a compensating  effect  in  the  lift  nonaveraging 
process  that  minimizes  the  influence  of  this  parameter  on  dispersion.  Other 
significant  results  of  the  analysis  are  as  follows: 

1.  Trajectory  deflection  occurs  in  a plane  that  leads  the 
plane  of  the  initial  angular  misalignment  by  approximately 
90  deg. 

2.  Dispersion  is  approximately  proportional  to  roll  rate  in 
contrast  to  lift  nonaveraging  dispersion  from  body-fixed 
configurational  asymmetries,  which  is  inversely  pro- 
portional to  roll  rate. 

3.  Dispersion  is  essentially  independent  of  vehicle  mass 
properties  with  the  exception  of  the  radius  of  gyration 
in  roll. 

4.  Dispersion  is  directly  proportional  to  the  angular 
misalignment  and  inversely  proportional  to  static  margin 
and  u sin  y. 
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APPENDIX  A. 


AVERAGE  VALUE  OF  THE  DISPERSION  INTEGRAL 


We  wish  to  find  the  average  value  of  the  integral 


rx 

= -iK / /x 


exp  f iK[,fn(  1 + x)  - x + 1 - fri2] } dx 


(A-l) 


as  X becomes  large.  Consider  first  the  real  part  of  .'/  and  make  the  change 
of  variable 


y = x - 1 


(A-2) 


such  that 


v = Re.?  = K / /rr  y sin  K |l  n (l  + \)  - y]dy  (A-3) 


We  define  the  average  value  of  v as  occurring  where  v"(y)  = 0,  as 
shown  in  Fig.  A-l.  From  Eq.  (A-3),  v,(y)  is  defined  by 


v'  (y)  = f(y)  sin  (y) 


(A-4) 


and  the  average  values  of  v(y)  occur  at  y = y where 


v"(yn)  = f'(yn)  si"  4Myn)  + f(yn)^  '(yn)  c°8  »Myn) 


= 0 


or 


%n> 

tan  ^ (yn)  = 

'n 


(A-5) 
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With  f(y)  = K y 1 + y)  and  4*  (y)  = ( 1 + ^ j - yj  from  Eq.  (A-  3), 

Eq.  (A-5)  becomes 

r / yn\  1 2K(l+y)2 

tan  il»  (yn)  = tan  K In  ( 1 + -j-;  - yn  = — (A-6) 

*-  J yn 

For  K sufficiently  large,  the  right  side  of  Eq.  (A-6)  is  positive  and 
» 1 such  that  the  angles  4<  (y  ) must  lie  in  the  first  and  fourth  quadrants 
near  ± it/2 


*iyn)  = - [(?°y '1)ir-  + pj  , n = 1,2,3,...  (A-7) 

where  p^  is  a small  angle.  Substitution  of  this  value  in  Eq.  (A-6)  and  the  use 
of  the  identity  for  the  tangent  of  a sum  yield 


tan  4<  (y  ) = x — - 
'7n  tan  i 


2K(1  + yn)‘ 
Z + y 


or,  because  p is  small, 
nn 


tan  p ^ p 
'n  ’ll 


2 + Yr 


2K(1  + yn)‘ 


( A-8) 


(A-9) 


The  values  of  the  precession  angle  4Myn)  where  the  average  values  of  v 
occur  are  then,  from  Eq.  (A-7), 


[ln(  1 

.1. 

(2n  - l)ir  , 2 + ^n 

y"J 

z ' z 

2K(l+yn)2 

(A-  10) 
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We  can  obtain  approximate  expressions  for  'Hyn)  an<3  f°r  Vn  by  using 
the  series  expansion  for  the  fn  term 


(A- 1 1 ) 


If  the  first  average  point  occurs  at  sufficiently  small  y^,  we  can 
approximate  £n(  1 + yj/2)  by  the  first  term  of  the  series.  With  this  approxima- 
tion and  neglecting  relative  to  n/2  for  sufficiently  large  K,  we  obtain  from 
Eq.  (A-  10) 


y j = */K 


(A-  12) 


The  net  transverse  velocity  from  Eq.  (A- 3)  with  the  first-order 
approximation  £n(l  + y/2)  « y/2  is 


v(Yl)  = K 


1 + y sin 


-2 


(A-  1 3) 


where  we  have  neglected  y relative  to  unity  in  the  square- root  term.  For 
subsequent  values  of  yn  sufficiently  small  that  the  approximations  that  con- 
tribute to  the  result  of  Eq.  (A- 13)  are  valid,  there  would  be  no  change  in  the 
subsequent  average  values  because 


Av 


yn^yn+l 


(2n+l)|  K 

2 f 2 sin  (^)  d (^)  = 0 

l2n-l)£ 


(A-  14) 
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It  can  be  shown  by  a procedure  similar  to  that  described  above  that 
the  imaginary  component  of  .^is  zero  for  all  sufficiently  large  K.  Hence,  the 
average  value  of  the  integral  ,^is 

- v + iw  = -2  (A-  15) 

which  is  independent  of  K for  sufficiently  large  K. 
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APPENDIX  B.  EFFECT  OF  DAMPING  ON  DISPERSION 


If  aerodynamic  damping  is  included  in  the  equations  of  motion,  Eqs.  (2) 
and  (3),  we  can  write  the  damping  moment  derivatives  M g and  M^  in  the  form 

= -v  (B- 1) 


where 


v 


T(Cm 


+ C ) 


m- 


(B-2) 


The  inclusion  of  damping  does  not  alter  the  quasi -steady  approximation 
for  precession  rates,  Eqs.  (8)  and  (11),  but  it  does  influence  the  angle-of- 
attack  convergence  in  the  form  (Ref.  3) 


— » (1  + 5)'  1 ^4exp(-b§)  = x’  1 ^2exp[-b(x2  - 1)]  (B-3) 


where 


pZIH 

r 


2u  mx  . sin  y 
st 


/ C + C ^ 
md 2 mq  mo’ 

1 - zi  v — c: 


'N 


( B-4) 


The  dispersion  integral,  Eq.  (22),  then  has  the  factor  exp[-b(x  - 1)] 
included  in  the  integrand 


= -iK  f Jx  exp[-b(x2  - 1)]  exp[i  (x)]  dx 


(B-5) 
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f 


f — r=* 


fS*JXDINO  MOB 


The  damping  is  so  small  at  the  high  altitudes  where  the  net  trajectory 
deflection  occurs  that  the  added  term  has  a negligible  effect  on  the  value 
of  the  integral  in  Eq.  (B-5).  This  becomes  apparent  if  we  change  the 
independent  variable  to  y defined  in  Eq.  (A-2),  which  gives 


//  = -iK 


L 


v/l  + y exp[-by(y  + 2)]  expfiiMy)]  dy 


(B-6) 


It  was  shown  in  Appendix  A that,  for  sufficiently  large  K,  the  first 

average  of  .if  occurs  at  y = Yj#  where  y^  « 1.  For  example,  from  the 

vehicle  and  trajectory  parameters  listed  in  Section  EL  with  C + C = -5, 

-4  mq  mo 

the  coefficient  b is  found  to  have  the  value  b = 2.05  X 10  . Thus,  the  damping 

term  exp[-by(y  + 2)]  remains  essentially  unity  over  the  integration  limits  and 

would  change  insignificantly  even  for  an  order  of  magnitude  increases  in  the 

damping  derivative  C 4-  C 

mm. 
q a 
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NOMENCLATURE 


aerodynamic  axial  force  coefficient 
aerodynamic  lift  force  derivative 
aerodynamic  normal  force  coefficient 
aerodynamic  normal  force  derivative 
cross-range  dispersion 
aerodynamic  reference  diameter 
altitude 
entry  altitude 

altitude  of  trajectory  deflection 

scale  height  for  exponential  atmosphere 

pitch  or  yaw  moment  of  inertia 

roll  moment  of  inertia 

radius  of  gyration  in  roll 

2Hpr/p  sin  y 

characteristic  length,  I/mxflt 

lift  force 

vehicle  mass 

pitch  moment 

pitch  damping  moment 

yaw  damping  moment 

roll  rate 

reduced  roll  rate,  pp/2 
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s 

t 

u 

u 

o 

V 

v,  w 

X 


z 


Y 

e 

e 

o 

0 

P 

5 

P 

Po 

Pi 

<t> 

4* 

*+ 

10 


aerodynamic  reference  area 
time 

vehicle  velocity 
entry  velocity 
dispersion  velocity,  v + iw 
components  of  V 

/TT1~ 


value  of  x at  trajectory  deflection 
static  margin 
altitude  variable 
path  angle 

angle  of  attack  (Euler  angle) 
entry  angle  of  attack 
pitch  rate 

V1 
2.  2 
<"  /Pr 

atmospheric  density 

reference  value  of  atmospheric  density 

atmospheric  density  at  trajectory  deflection  altitude 

roll  angle  (Euler  angle) 

precession  angle  (Euler  angle) 

precession  rate 

precession  modes 

undamped  natural  pitch  frequency 


4 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high -power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro- optic  a;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmoephere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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